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SYSTEM AND METHOD FOR DETERMINING CAPACITANCE 
FOR LARGE-SCALE INTEGRATED CIRCUITS 



TECHNICAL FIELD OF THE INVENTION 

The present invention is directed, in general, to integrated 
circuit (IC) design and, more specifically, to a system and method 
for determining capacitance for large-scale ICs. 

BACKGROUND OF THE INVENTION 

Interconnect analysis for modern ICs has traditionally 
required large scale capacitance calculations. However, such 
analysis required the geometry of the ICs to be described in 
detail, which resulted in millions, or even billions, of polygonal 
elements. The interconnect analysis yielded a sparse capacitance 
matrix of dominant coupling coefficients between conductors 
(commonly called "nets") . A net is formed of metallic pieces that 
are located within the IC and coupled together such that they 
cooperate to conduct current between two points. The net also has 
a capacitance associated with it. 

Traditionally, capacitance calculations have not required 
great accuracy (accuracies of 10-20% were acceptable) . To achieve 



this relatively low accuracy, pattern matching schemes consisting 
of precharacterized libraries were used. However, nets such as 
clock nets, power grids, sensitive analog nodes or critical paths 
of a circuit require a higher accuracy in their solution. 
Unfortunately, as the requirement for accuracy increased, pattern 
matching solutions became insufficient. 

Numerical approaches were then explored in an effort to 
increase the accuracy in capacitive calculations. Conventional 
numerical approaches rendered conductor surfaces (nets) in terms of 
discrete elements. Enforcing the equation either as a set of 
collocation points or as a distribution, as with a Galerkin scheme, 
led to a dense linear system. The large size of these systems 
usually required that they be solved with Krylov iterative methods. 
The required matrix-vector products were accelerated with 
compressed-matrix schemes (such as FastCap) # pre-corrected Fast 
Fourier Transforms, or Singular Value Decomposition ( "SVD" ) -based 
representations. While these methods are capable of high accuracy, 
they are restricted to systems of only a few hundred thousand 
unknowns. Even with only one unknown per polygon, large problems 
having millions of polygons remained well beyond the reach of these 
methods . 

The industry further tried what is conventionally termed a 
"random walk" method, based on the well-known Monte Carlo method. 



Unfortunately, the random walk method required a significant amount 
of central processing unit ("CPU") time to achieve high accuracy. 
For this reason, larger problems involving millions of polygons, 
proved to be unsatisfactory. 

Next, finite-difference or finite-element methods were 
attempted. These methods required that a three-dimensional mesh be 
imposed on the geometry of the capacitive structure. However, 
these methods created more unknowns than did other numerical 
methods, and therefore were even less desirable. 

Accordingly, what is needed in the art is way to increase the 
accuracy of IC capacitance calculations without incurring excessive 
computational time or difficulty. 




SUMMARY OF THE INVENTION 

^ *$p address the above-discussed deficiencies of the prior art, 
the prestent invention provides, for use in an integral equation 
f ormulatiori of capacitance, a system for, and method of, generating 
a representation of charge distribution for a given capacitive 
structure (whicfi may be an integrated circuit) . In one embodiment, 
the system included: (1) a charge variation function generator that 
creates a multidimensional charge variation function that is 
independent of a conductive geometry of the structure and (2) a 
conductive geometry generator, associated with the charge variation 
generator, that creates a conductive geometry that is independent 
of charge variation in the structure, the charge variation function 
and the conductive geometry employable in the integral equation 
formulation to reduce a complexity thereof . 

The present invention therefore^ introduces the broad concept 
of rendering the charge variation \and conductive geometry 
independent of one another, thereby significantly reducing the 
complexity of the subsequent integral equation^f ormulation required 
to determine a structure's capacitance. 

In one embodiment of the present invention, the integral 
equation formulation is a Fast Distribution Method. Those skilled 
in the pertinent art will understand, however, that the broad scope 



of the present invention is not limited to a particular integral 
equation formulation. 

In one embodiment of the present invention, the charge 
variation function is a three-dimensional function. Three- 
dimensional functions are particularly useful in describing 
integrated circuit structures. However, the present invention is 
not limited to a particular dimensionality. 

In one embodiment of the present invention, the charge 
variation function is a smooth function of spatial location. 
Smooth functions can be accurately represented by low-order 
polynomials and are therefore advantageous. 

In an illustrated embodiment of the present invention, the 
conductive geometry generator iteratively creates the conductive 
geometry. In the illustrated embodiment of the present invention, 
the charge variation function generator employs a generalized 
minimal residual -based krylov method to determine the 
multidimensional charge variation function. 

In one embodiment of the present invention, the conductive 
geometry is represented in an octtree. The structure and function 
of octtrees will be set forth in the Detailed Description that 
follows. Those skilled in the pertinent art will understand, 
however, that octtrees are not required by the present invention. 

The foregoing has outlined, rather broadly, preferred and 




alternative features of the present invention so that those skilled 
in the art may better understand the detailed description of the 
invention that follows. Additional features of the invention will 
be described hereinafter that form the subject of the claims of the 
5 invention. Those skilled in the art should appreciate that they 
can readily use the disclosed conception and specific embodiment as 
a basis for designing or modifying other structures for carrying 
out the same purposes of the present invention. Those skilled in 
!S » the art should also realize that such equivalent constructions do 
if not depart from the spirit and scope of the invention in its 
: 4 broadest form. 

} _£ 

: 
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BRIEF DESCRIPTION OF THE DRAWINGS 



For a more complete understanding of the present invention, 
reference is now made to the following descriptions taken in 
conjunction with the accompanying drawings, in which: 

FIGURE 1 illustrates a capacitive structure that forms one 
environment with which the present invention can operate on to 
reduce the complexity thereof; 

FIGURE 2 illustrates a block diagram a systetnfor__det^rmining — • 
a capacitance of an djntegxa^redT^Trcuit constructed according to the 
pw^TcipTes of the present invention; 

FIGURE 3A illustrates a method of refining the approximation 
of the geometry of a capacitive structure using octtrees; 

FIGURE 3B illustrates the refined octtree representation of 
the geometry of a capacitive structure of FIGURE 3A; and 

FIGURE 4 illustrates a table of the time and memory required 
to determine the capacitance of a section of a digital chip 
employing an embodiment of an algorithm of the present invention. 
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DETAILED DESCRIPTION 

Referring initially to FIGURE 1, illustrated is a capacitive 
structure that forms one environment with which the present 
invention can operate to reduce the complexity thereof. The 
5 capacitive structure 100 illustrates an example of three nets CI, 
C2, C3 in an integrated circuit (IC) . A net is comprised of 
metallic pieces connected together to allow a current to pass from 

:~l a first point in an IC to a second point in the IC. The net has an 
associated capacitance. The net may pass through and include 

IS; portions of different IC structures that can affect the associated 
capacitance of the net. The IC structures may include traces, 
capacitors, switches, resistors and other elements that make up an 

1 IC " 

^ The net CI is a block representation of a first net within the 

15 capacitive structure 100. The net C2 is a block representation of 
a second net within the capacitive structure 100. The net C2 is 
associated with, but not coupled to, the net CI. The net C3 is a 
block representation of a third path within the capacitive 
structure 100 and is associated with, but not coupled to either one 
20 of, the net CI and the net C2 . 

The capacitive value of the capacitive structure 100 depends 
upon the geometry of each of the nets CI, C2, C3, the capacitance 
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for each of the nets CI, C2, C3 and the electromagnetic interplay 
between the nets CI, C2, C3 . The capacitive value also depends 
upon the charge distribution for the capacitive structure 100. 

As previously described, the difficulty with previous 
capacitance calculation approaches is that they require the 
geometry of the nets to be captured accurately, causing the minimum 
problem size for the capacitance calculation to exceed practical 
limits. The present invention addresses this substantial 
shortcoming by providing a method that does not require that the 
geometry of the nets be accurately captured. 

Consider the capacitance of the net CI; it is affected by the 
net C3 . However, the net C3 is located far enough away that the 
exact geometry of the net C3 is not critical. Therefore, were the 
net C3 to be replaced with a simpler geometric description, the 
number of unknowns could be reduced. In particular, the number of 
unknowns needed to determine the capacitance of the capacitive 
structure 100 would then be constrained by the charge distribution, 
rather than the geometry, of the nets. 

In one embodiment of the present invention, the charge 
variation is decoupled from the geometry of the nets. Decoupled 
means the charge variation is determined separately from the 
geometry. However, the charge variation and the geometry are used 
in the actual capacitance calculation. Also, the surface charge 
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distribution is represented as a smooth function of spatial 
location which is projected onto the geometry. This results in the 
number of unknowns being dictated solely by the charge variation. 

The geometry of the nets CI, C2 , C3 is described 
hrer^xchically and captured to a level of detail needed to 
determine an accurate^golution. For example, in FIGURE 1, the net 
CI and the net C2 are captured with rela^El/v^ty-^det^ailed geometric 
descriptions. However, the net C3 is captured with relatively- 
coarse geometric descriptions. 

With the charge variation decoupled from the geometry of the 
nets CI, C2, C3, the capacitance is calculated using an iterative 
linear solution. In one embodiment of the present invention, the 
present invention uses a Krylov method based on a Generalized 
Minimal Residual (GMRES) method to solve the capacitance of the 
nets CI, C2, C3 . For more background material on this subject, 
see GMRES: A Generalized Minimal Residual Algorithm for Solving 
Nonsymmetric Linear Systems , by Y. Saad and M. H. Shultz, SIAM 
Journal on Scientific and Statistical Computing, 7(3) : 856-69, 1986 
(incorporated herein by reference) . 

Unlike standard GMRES, the present invention represents the 
net geometry as a set of the charge distribution functions rather 
than a set of discrete vectors. This method of representation is 
founded on the assumption that the charge variation function is a 




smooth function of spatial location. With the advantageous use of 
the charge distribution functions, the present invention makes it 
possible to solve problems that are intractable with conventional 
methods . 

Turning now to FIGURE 2, illustrated is a block diagram a 
systeoifor determining a capacitance of an IC, generally designated 
200, constructed— according to the principles of the present 
invention. The system 200^^eoin^rises an integral equation 
formulator 210, a charge variation function>g^nerator 220, and a 
conductive geometry generator 230. The integr^rl^equation 
formulator 210 determines the total capacitance of all the net^C^, 
C2, C3 within the IC. 

In determining the capacitance for the nets CI, C2, C3 , the 
integral equation formulator 210 uses an iterative linear solution 
method. The iterative linear solution method is a recursive method 
that starts with a first guess of the geometry of a net and 
converges toward a charge distribution function that is within 
acceptable limits when applied to the geometry. If the charge 
distribution function is not within acceptable limits, the geometry 
is divided into smaller geographic elements called subdivisions. 
The iterative linear solution method is then performed on the 
subdivided structure. The resolution of the charge distribution 
function and the subdivision of the geometry continues until the 
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charge distribution function for each of the subdivisions is within 
acceptable limits. 

Once the integral equation formulator 210 has found the 
solutions for the charge distribution functions for each of the 
nets CI, C2, C3, the integral equation formulator 210 determines 
the capacitance for the nets CI, C2, C3 . 

Associated with the integral equation formulator 210 is the 
charge variation function generator 220. The charge variation 
function generator 220 creates charge variation functions that the 
integral equation formulator 210 uses to refine the representation 
of the charge distribution. For each iteration of linear solution 
method, the charge variation function generator 220 creates a 
charge variation function for a particular subdivision of the 
geometry. 

Also associated with the integral equation formulator 210 is 
the conductive geometry generator 230. The conductive geometry 
generator 230 generates an advantageous description of the geometry 
of the IC by performing an adaptive subdivision of the nets CI, C2, 
C3 . Over the course of the iterative linear solution method, the 
conductive geometry generator 23 0 may further subdivide the 
existing subdivisions of the geometry to assist the integral 
equation formulator 210 in reaching an acceptable solution. 

In ^The- — illuptra £e d _^embodiment , the integral equation 
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^ormulator 210 uses a two -component representation of the charge 
distribution. The first component, x R / is the characteristic 
function of^he surfaces of the net and represents the geometry, 
without regard to the charge distribution. The second component is 
a charge distribution ^function f that represents the charge 
distribution, without regard tb-<he geometrical layout. Since the 
charge distribution and the geometry\are wholly decoupled, the 
representation of the charge distribution prbg^cted on the geometry 
requires far fewer unknowns than previous method^ 

In previous numerical analysis methods, the geometry of the 
entire set of net surfaces were incorporated into the calculations 
for any subdivision of the geometry. As a result, these methods 
required solving for a relatively large number of unknowns to 
account for all of the geometry of the nets CI, C2, C3 . This 
included all the geometry that was far-removed from the subdivision 
being considered. In contrast, the present invention 

advantageously does not require the incorporation of all the far- 
removed subdivisions. This greatly reduces the number of unknowns 
required in solving for capacitance of the IC. 

In one embodiment of the present invention, the characteristic 
geometry func^orT^^^JUieset of net surfaces, x R / is a predefined 
function. In an alternate embodimentoT^he^present invention, x R 
is calculated from information related to the capacitive sFrtlcfctuce 
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of the\ IC. Xr is scalar-valued function that has a n-dimensional 
positions, as the input. To illustrate the form of x R / suppose the 
geometry tb be described is the unit square in the xy-plane. The 
form of Xr foqr this geometry then would be: 
\ Zr = s(x)s(y)S{z) 9 

where s (x) =1 for O^x^l and s(x) is zero elsewhere, and 5(z) is the 
standard Dirac delfca function. Those skilled in the art are 
familiar with the properties and use of Dirac delta functions. 

The function f, wrtiich represents the charge distribution 
without regard to geomeory, is evaluated over R, the set of 
surfaces of the nets CI, C2,\c3. The surface integrals over R are 
reformulated as volume integrals via the relation: 

lf\Lfz* 

where V is the volume that the net Comprises. 

To solve the iterative linear solution, the integral equation 
formulator 210 uses an initial guess g, for the charge 
distribution, and a starting subdivision of the geometry of the 
net. The integral equation formulator \210 also creates the 
projection matrix P R from x R , \ 

k 



-14- 




\ 

The\^bracketed expression above is the standard inner product for 
functions : 



Also in this equation, P k is a primitive projection matrix where 
the (i,j) entry of P k = ^liljl k . Finally the liS represent the choice 
of orthogonal polVnomials. In one embodiment of the present 
invention, LegendreXpolynomials form the system of orthogonal 
polynomials. One skilled in the pertinent art is familiar with and 
the use of Legendre polynomials. Also, P b R will represent that part 
of P R that contains information of the geometry of subdivision b. 

Using the initial guess gNCor the charge distribution and the 
starting subdivision of the geometry of the net, the integral 
equation formulator 210 determines charge distributions for a given 
potential on the nets CI, C2, C3 . Frbm the charge distributions, 
the integral equation formulator 210 determines the capacitance of 
the nets CI, C2, C3 . 

In the determination of the charge distributions, the 
iterative linear solution uses the following fundamental equation 
that relates the charge and the potential: 
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wh^re G (r, r' ) is the Green's function for the geometry of the IC 
and p\r' ) is the charge density. Those skilled in the art are 
f amiliarNvith the properties and use of Green's functions. 

In one\embodiment of the present invention, the fundamental 
equation relating the charge and the potential is expressed as: 

W= MP R f 

where the \|i is a potential distribution and f represents the charge 
distribution without regard to the geometry. From this equation, 
the integral equation fo^mulator 210 can determine the charge 
distribution function f. 
v In one embodiment of t&e present invention, the integral 
equation formulator 210 uses a KWlov method based on a variation 
of a Generalized Minimal Residual Method (GMRES) for the iterative 
linear solution. Background information concerning GMRES is 
discussed in GMRES: A Generalized Minimal Residual Algorithm for 
Solving Nonsymmetric Linear Systems , by y\ Saad and M. H. Shultz, 
SIAM Journal on Scientific and Statistical Computing, 7(3) : 856-69, 
1986 (incorporated herein by reference) . 

GMRES is used to solve for the charge distribution function f 
in the following equation: 
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or, equivalently : 



P R MP R f=P R f Q , 



where % is a\potential which is 1 over an individual net and 0 
elsewhere. Tho^se skilled in the pertinent are familiar with the 
application of GMRES . Nevertheless, below is a description of the 
process both to give a complete picture of the iterative linear 
solution process and to show the advantages of decoupling the 
charge variation fromVthe geometry. 

The integral equation formulator 210 starts iterative linear 
solution with an initial ^uess g for the charge distribution and an 
initial geometry. The integral equation formulator 210 also 
determines x b R / which represents the geometry of the box b in the 
subdivision, and P b R , which isNthe part of the projection matrix P R 
which covers the box b (see\FIGURE 3A for more information 
concerning geometry subdivisions and boxes) . From this 

information, the integral equation formulator 210 computes a 
potential i|r based upon the initial change distribution guess g and 
the current subdivision, 

The integral equation formulator ^10 also determines the 
desired potential \Ji 0 : 

where o is the constant function 1 and P s is defined as P R is over 
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an individual net S. The integral equation formulator 210 then 



determines the difference between the desired potential ty 0 and the 
potential's This difference, ty Q -ty t becomes first charge variation 
function f x . \Note that f x is not a representation of the charge 
distribution described by f . ^ is a function used to modify the 
charge distribution function f to reach a desired resolution. 

The integral equation formulator 210 determines the ratio y/3/ 
where y = HfjR and p =\lh|f 0 llR. If the ratio is sufficiently small, 
the iterative linear solution is complete. If the ratio is 
unacceptably large, the intsegral equation formulator 210 normalizes 
f 1 and proceeds with the linear iterative solution. 



V The integral equation formulator 210 then invokes the charge 
variation function generator 220 to create a charge variation 
function which refines the description of the charge distribution. 
At the beginning of the iterative\process the charge variation 
function generator 220 creates a second charge variation function 
f 2 . In general, the formula for the crfeation of the j+l th charge 
variation function from the j th charge variation function is: 



\ 




= MPifj . 



Where M is an operator that converts a charge density to a 
potential distribution. The (i,j) entry \of the matrix 

representation of M is : \ 



Where l b ± \s the i th moment in box b. Note that the operation of M 
involves al\ boxes, not just box b. G(r,r') is the Green's 
function for the set of net surfaces R. Those skilled in the art 
are familiar wioh the properties and use of Green's functions. 

The charge variation function generator performs the operation 
M using a variant V>f the Fast Multipole Method (FMM) algorithm 
called the Fast Distributed Method (FDM) . The FDM differs from the 
FMM in several areas. Virst, the FDM omits the FMM' s initial step 
of computing multipole representations from point charges. FDM 
omits this step because theV input is already in terms of a charge 
distribution. Second, the FDM omits FMM' s direct point-to-point 
interactions. Instead, the FDMVuses a source charge distribution 
in a cube. The potential distributions are calculated in all 
neighboring cubes, including the source cube itself. 

Third, the FDM omits FMM's final\step of interpolating point 
potentials from the local expansions. Tne FDM omits the last step, 
because the desired output is a potential ^di str i^ut ion. Finally, 
FDM uses Legendre polynomial expansions for both charge 
distributions and local expansions. The \FMM uses multipole 
expansions for representing charge distributions and local 
expansions for representing potential distributions. 
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The FMM algorithm, Legendre polynomial expansions and their 
use arev well known to those skilled in the pertinent art. 
Background\ information concerning numerical analysis and 
capacitance calculations is discussed in Introduction to Numerical 
Analysis . by J. \Soer and R. Bulirsh, Springer-Verlag 1979 and in 
Preconditioned, Adaptive, Multipole-Accelerated Iterative Methods 
for Three -Dimensional First-Kind Integral Equations of Potential 
Theory , by K. NaborsX et al., SIAM Journal on Scientific and 
Statistical Computing, 1S\(3) : 713-735, May 1994 (both incorporated 
herein by reference) 

The use of the FDM results in a significant savings in time 
over both the analytical calculation using the explicit form of M 
and the FMM. The savings are dueVto the omission of several FMM 
steps and the omission of FMM' s direc\ point-to-point interactions. 

Once the charge variation function generator 220 has created 
the charge variation function f j+1 , the integral equation formulator 
210 uses this new function to create th^e (i, j) entry of a 
Hessenberg matrix: 



H ij-{fi>fj+i) } 



for all i< j . 

The integral equation formulator 210 uses the Hessenberg matrix to 
solve for the new coefficients of the charge variation ^functions to 
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be used \to further refine the charge distribution function f. 

Af teis the entries in the Hessenberg matrix are computed, the 
integral equation formulator 210 then orthogonal izes f j+1 from the 
other fiS. *Ehe charge variation function generator 220 then 
creates the (j\l,j) entry of the Hessenberg matrix by the 



expression: 



H 



f, 



This expression is the Nstandard norm of functions defined by: 



After the creation of this last entry of the Hessenberg matrix, the 
integral equation formulator 2^0 then normalizes f j+1 . The 
orthogonalization and the normalization of f j+1 create orthonormal 
charge variation functions. As a Result, the integral equation 
formulator 210 is able to obtain the Correct charge distribution 
function f very quickly. 

The integral equation formulator 21(3 then solves for x, the 
coefficient vector for the charge variations, using least-squares 
on the equation: 

Hx - ye x 



-21- 



• # 

where H is the Hessenberg matrix generated by the integral equation 
formulator 210 as explained above. e x is the 1 st unit vector with 
the forta (1,0,0,0,0, . . .) (with j number of zeros) . The vector x 
is the coefficient vector for the charge distributions and y is the 
norm of theXfirst charge variation function, y is defined as: 

\ r = ll/.IL • 

Since the charge Variation and the geometry are decoupled, the 
representation of the* f t will not become very large as compared to 
the previous methods. tJsing the residual of the least squares, the 
integral equation formulator 210 computes the factor r as the 2- 
norm of the least squares residual. 

The integral equation rormulator 210 then generates a new 
guess for the charge distribution function f using the equation: 

i=\\ 

The integral equation formulator 2r0 then determines the ratio 
r/p. Where r is the 2-norm of the least Nsquares residual r and (3 
is the norm of the desired potential (defined, earlier as (3 = lk|/ 0 llR) . 
If the ratio is within acceptable limits, the^charge distribution 
has converged and the integral equation formulator 210 terminates 
the iterative linear solution. \ 
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If\the ratio is not within acceptable limits, the integral 
equation \formulat or 210 uses the new guess for the charge 
distribution function f, calculated above, as the starting point 
for the next iteration of the iterative linear solution. However, 
before the integral equation formulator 210 performs the next 
iteration, the integral equation formulator 210 must determine if 
the current subdivision requires a more refined approximation. 

The integral equation formulator 210 determines if the current 
subdivision b requires^ refinement by calculating a charge-geometry 
error e. In one embodiment of the present invention, the integral 
equation formulator calculates the charge -geometry error e b 
associated with box b in \the subdivision using the following 
equation: 

,N~ 1 



Where li is the i th moment basis function. f represents the 
distribution of the charge without any reference to the geometry of 
the IC and f ± b is the component of f alonisj the i th moment for the 
charge function contained in the subdivision b. The summation of 
1 ± between k and N-l represent the high\ order polynomials. 
Finally, N is the number of moments used in thfe expansion of f . 

The integral equation formulator 210 then determines if the 
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charge -geometry error is within acceptable limits using a charge- 
jeometry error criterion. The charge -geometry error criterion is 
defined as : 



R 



Where S is the minimal subdivision containing the selected net. e 
is a small constant that determines the allowable tolerance of the 
criterion. Also, ||e b || R is the magnitude of the charge -geometry 
error in subdivision b. Note that the ers^or threshold is relative 
to the charge distribution f only on they subdivision of the 
selected net. 




If the charge -geometry error criterion is within acceptable 
limits, the integral equation formulator 210 has finished the 
resolution of the charge variation function. If the charge- 
geometry error criterion is not within acceptable limits, then the 
integral equation formulator 210 calls the conductive geometry 
generator 230 to refine the approximation of the geometry by 
subdividing the current subdivision b. See FIGURE 3 for a 
discussion of the refining the approximation of the geometry of a 



net . 



After the conductive geometry generator 23 0 has subdivided the 
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geometry, the integral equation formulator 210 then starts the 
iterative linear solution again on the new subdivision. The 
integral equation formulator 210 uses the last calculated charge 
distribution function f as the starting point for the calculations 
on the new subdivision. 

Once the integral equation formulator 210 has completed the 
iterative linear solution calculations for all the nets, the 
integral equation formulator 210 then calculates the capacitance of 
the IC. To determine the capacitance, the total charge on each net 
is computed from the charge distribution. The charge on the net 



to a single method of determining capacitance. Other methods of 
determining capacitance from charges on nets are within the broad 
scope of this invention. Background information concerning the 
calculation of the capacitance is discussed in Classical 
Electrodynamics , by John Jackson, John Wiley & Sons (3d edition, 
1998) (incorporated herein by reference) . 

Turning now to FIGURE 3A, illustrated is a method of refining 
the approximation of the geometry of a capacitive structure (or 
net) using octtrees. An octtree is a generalized way of 
subdividing the geometry into rectangular "boxes" that are 
congruent within every level of subdivision. Each box in a 
subdivision is called a "leaf." 



whose surface is S is 




is 



The present invention is not limited 




The conductive geometry generator 23 0 starts with the initial 
geometry 310 of a net. When the integral equation formulator 210 
determines that the initial geometry 310 needs to be refined, the 
conductive geometry generator 230 divides the initial geometry 310 
5 into three leaves shown in the first subdivision 320. Usually the 
geometry generator 23 0 will divide a given portion of a geometry 
into four equal leaves. However, the conductive geometry generator 
230 will discard any leaf that is not required for the 
^ approximation of the geometry. In this instance, the conductive 
iM geometry generator 230 discarded the lower right leaf in first 
LI subdivision 320. 

ijjj On the next iteration of refining the approximation of the 

geometry, the conductive geometry generator 230 subdivides the 
U lower left leaf, of the first subdivision 32 0, into three smaller 
16 leaves shown in the second subdivision 330. In this iteration, the 
^ r conductive geometry generator 230 discarded the upper right leaf of 
the second subdivision 330. 

On the final iteration of refining the approximation, the 
conductive geometry generator 230 subdivides the upper right leaf, 
20 of the first subdivision 320, into three smaller leaves shown in 
the third subdivision 340. In this iteration, the conductive 
geometry generator 23 0 discarded the lower right leaf of the third 
subdivision. 
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Now turning to FIGURE 3B, illustrated is the refined octtree 
approximation of the geometry of the capacitive structure of FIGURE 
3A. The refined approximation 350 shows the results from the 
conductive geometry generator subdividing the initial conductive 
geometry 310. Note that the leaves are not of equal size. In 
particular, the leaf on the upper- left quadrant of 3 50 is four 
times larger than the other leaves illustrated. This is a result 
of the fact that although the n-dimensional boxes at a particular 
level are all the same size. However, all the leaves are not 
required to be the same size since the leaves are not generated 
from the same level of subdivision. The subdivision of one leaf 
does not mandate the subdivision on all other leaves that are the 
same size. 

There are various advantageous reasons for using octtrees to 
approximate the geometry of a conductive structure. First, the 
description of any n-dimensional space is relatively simple. 
Second, octtrees provide a simple description of the geometry. 
Third, the boxes (leaves) that comprise an octtree fit together in 
a natural way regardless of how often they are subdivided. Fourth, 
it is a relatively simple task to ignore sections that are not 
important to the approximation of the geometry. Finally, any 
octtree can be refined to an arbitrary level of accuracy. 

Turning now to FIGURE 4, illustrated is a table 400 of the 



time and memory required to determine the capacitance of a section 
of a digital chip employing an embodiment of an algorithm of the 
present invention. The algorithm was performed on an ULTRASparc II 
workstation (manufactured by Sun Microsystems) having a 400 
Megahertz (MHZ) CPU and 1 Gigabyte of memory, which is a typical 
system that may be employed in performing the capacitance 
calculations of the present invention. The algorithm used octtrees 
with square subdivisions as the geometrical structure and employed 
a Krylov algorithm to solve for the charge variation. The 
tolerance of the system was set at about 5%. 

The time and memory requirements illustrated are derived from 
calculations of an interconnect for the section of the digital 
chip. The section is about 0.5 mm on a side and is fabricated 
using 0.25 micron complementary metal oxide semiconductor ("CMOS") 
technology. The section contains 3,400 nets and geometry is 
described with 600,000 rectangles. Solving for the capacitance of 
a critical net required approximately two hours of CPU time and 610 
Megabytes of memory. The calculation required a large amount of 
memory, in part, because this critical net of the digital chip is 
very close to many other nets. 

The table 400 illustrates an exemplary time required by each 
part of the algorithm at various levels of refinement. The table 
400 includes columns for (1) the number of cubes, (2) the time per 
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iteration for the Fast Distribution Method, (3) the geometry 
manipulation time, including miscellaneous operations, (4) the 
total time per iteration and (5) the memory used. At the finest 
level there were approximately 5.8 million unknowns (with 35 
moments for each cube) . The algorithm required about 100 seconds 
to perform a FDM calculation. 

The present invention thus allows companies to solve problems 
many times larger than possible using other conventional methods. 
The companies may therefore perform computational tests of proposed 
ICs without having to construct the proposed IC. The present 
invention provides the end user with the benefits of flexibility 
and speed. Flexibility may result from a choice of methods for 
solving for both the charge variation and the advantageous 
geometry. Speed is provided by the categorical division of the 
problem between solving for the charge variation and solving for 
the advantageous geometry. 

Although the present invention has been described in detail, 
those skilled in the art should understand that they can make 
various changes, substitutions and alterations herein without 
departing from the spirit and scope of the invention in its 
broadest form. 



